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Abstract Electrostatic potential of electrons near the surface of static strange 
stars at zero temperature is studied within the frame of the MIT bag model. We 
find that for QCD parameters within rather wide ranges, and if the nuclear crust on 
the strange star is at a density leading to neutron drip, the electrostatic potential 
is insufficient to establish an outwardly directed electric field, which is crucial for 
the survival of such crusts. If a minimum gap width of 200 fm is called in for a 
more stringent constraint, our calculations completely rule out the possibility of 
such crusts' presence on strange stars. Therefore, our results prefer against the 
existence of neutron-drip crust in nature. 

Key words: dense matter - elementary particles - stars: neutron 
1 INTRODUCTION 

There may exist a more stable state of hadrons than 56 Fe, called strange quark matter (SQM), 
which is a bulk quark phase consisting of roughly equal numbers of u, d, and s quarks plus 
a smaller number of electrons to guarantee charge neutrality (Witten 1984). This theoretical 
result is unlikely to be proved or disproved through QCD calculations at least in the foreseeable 
future. Final adjudication must come from experiments conducted with accelerators or from 
astrophysical tests. One of the important consequences of Witten's hypothesis is the prediction 
of strange stars (Haensel, Zdunik, & Schaeffer 1986; Alcock, Farhi, & Olinto 1986, hereafter 
AFO), i.e., stars made of SQM. The presence of such stars should never be rare, in fact some 
authors even inferred that all neutron stars might have been converted to strange stars, since the 
whole Galaxy is likely to be contaminated by stranglets (Glendenning 1990; Madsen & Olesen 
1991; Caldwell & Friedman 1991; Medina- Tanco & Horvath 1996). For decades, unremitted 
efforts have been made to observationally discriminate strange stars from neutron stars, and 
some candidates for strange stars have been accumulating (Cheng et al. 1998; Li et al. 1999; 
Xu et al. 2001a, b). However, even to now it is still premature to reach any firm conclusions. 
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A strange star can have a normal matter crust. The density at the base of the crust is a 
crucial parameter, which has been discussed by many authors (AFO; Huang & Lu 1997a, b; 
Xu & Qiao 1999; Yuan & Zhang 1999; Chen & Zhang 2001; Ma et al. 2002). Recently, Ma et 
al. (2002) studied the influence of some parameters, such as the strange quark mass (m s ), the 
bag constant (B), and the strong coupling constant (oe c ), on the crust density. However, their 
study does not give out clearly the ranges of parameters that enable a crust as dense as neutron 
drip density. Here we study the problem in more detail. The parameter ranges that support 
a neutron-drip-density crust will be calculated numerically. Our paper is arranged as follows. 
Section 2 describes the existence of strange star crust. Section 3 gives the detailed procedure 
of our calculations and the numerical results in the a c — m s — B space. Finally, we summarize 
the results in Section 4. 



2 THE EXISTENCE OF CRUST 

It was pointed out by AFO that a strange star could be covered by a normal material crust. The 
presence of electrons in SQM is vital to the existence of such crusts. Because s quark's mass 
is larger than u and d quarks, it is slightly deficient in equilibrium SQM. Electrons are thus 
called in to make the system electrically charge neutral. As quarks are bounded through strong 
interaction, they should have a very sharp surface with thickness of the order of 1 fm. On the 
other hand, the electrons, bounded by the Coulomb force, can extend several hundred fcrmis 
beyond the quark surface. So, a strong electric field, ~ 10 17 V cm -1 and outwardly directed, 
will be established in a thin layer of several hundred fermis thick above the strange matter 
surface. This field can support a crust, composed of normal nuclear matter, suspended out of 
contact with the SQM core. 

The so-called neutron drip density, pdrip ~ 4.3 x 10 11 g cm~ 3 sets an absolute upper limit 
for the density of the crust, p cr ust- The reason is that, if p crU st reaches pdrip, neutrons will 
begin to drip out. Being electrically neutral, the neutrons will fall freely into the core, and, by 
hypothesis, be deconfmed to be SQM. As a consequence, p crus t will keep going down, until it's 
below pdrip- Conventionally, when the nuclear crust is taken into account, the bottom density 
was assumed to be pdrip- However, is neutron drip the only limit on the crust density? Could 
the crust density actually go so far? Huang & Lu (1997a, b) said NO to both questions. By 
proposing that mechanical balance should be held between electric and gravitational forces 
on the whole crust, not only on a single nucleus as modelled by former authors (AFO), they 
claimed that at a density still far lower (~ Pdrip/5) than the neutron drip density, the crust 
would begin to break down. 

AFO proposed a model to describe the gap between the SQM core and the nuclear crust 
within the framework of Thomas-Fermi model. The electric field should be described with the 
following Poisson's equation: 

' 4a(V 3 - V 3 )/[3ir(hc) 2 }, z < 0, 
= J 4aV 3 /[M^)% 0<z<z G , (1) 

k 4a{V 3 - V 3 )/[3ir(hc) 2 }, z G < z, 

where z is a space coordinate measuring the height above the quark surface, a is the fine- 
structure constant, V 3 /Sir 2 h 3 is the quark charge density inside the quark matter, V is the 
electrostatic potential of electrons, V c is the electron Fermi momentum deep in the crust, which 
represents the positive charge density of the ions within the crust, and zq is the gap width 
between the SQM surface and the base of the crust. In fact, V q and V c are the boundary values 
of the above equation: V — > as z — > — oo, V — > V c as z — > +oo. Meaningful solutions to the 
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above equations occur only ii V c < V q . Actually in a core-crust system, V q turns out to be the 
electron electrostatic potential at R m . Here R m represents the maximum radius below which 
electrical charge neutrality is locally satisfied in the SQM core. Certainly R m is smaller than 
R, the radius of the core. 

Since the electron chemical potential (we will prove that it equals to V later) at which 
neutron drip occurs is ^26 MeV (Baym, Pethick, & Sutherland 1971), to keep a crust at a 
density of pdrip suspended, the electrostatic potential of electrons near the edge of the SQM 
core, V q , must at least be larger than 26 MeV. In this work, we calculated V q for static SQM 
cores at zero temperature. Because of the uncertainties inherent in the critical QCD-related 
parameters, the s quark mass, m s , the strong interaction coupling constant, a c , and the bag 
constant, B, our calculations were actually carried out in a space expanded by those three 
parameters. We found that for parameters within rather wide ranges, to support a crust at a 
density leading to neutron drip is not possible. Especially, for the conventional choice of the 
parameters, i.e., m s =200 MeV, a c =0.3, and SV 4 =145 MeV, our calculations indicate V q ~20 
MeV, well below 26 MeV. 

It is interesting to mention that the properties of charm-quark stars (viz. strange-quark 
stars with an additional charm-quark population) have been studied by Kettner et al. (1995). 
But unfortunately, they found charm-quark stars are unstable against radial oscillations, i.e., 
no such stars can exist in nature. So our work will only pivot around strange-quark stars. 



3 ELECTROSTATIC POTENTIAL OF ELECTRONS IN SQM CORES 

3.1 Governing Equations 

The property of SQM is generally described using the phenomenological MIT bag model 
(Chodos ct al. 1974), which simplifies the dynamics of confinement by introducing an approxi- 
mation that the quarks are separated from the vacuum by a phase boundary and the region in 
which quarks live is endowed with a constant universal energy density B. Since the description 
of SQM has been introduced elsewhere (Farhi & Jaffe 1984; Kettner et al 1995), we will go 
straightforward to the governing equations and describe the procedure of our calculations. Our 
goal is to determine the electrostatic potential of electrons at R mi i.e., V q . 

We assume that charge neutrality is locally satisfied. This is true popularly: Due to the 
rearrangement of electron charge inside and outside of the surface of an SQM core, the redundant 
positive charge of the quarks will be balanced locally by e~ only up to radial distances r < R m . 
However, R m in fact is only minutely smaller (several hundred fermis, see AFO; Kettner et al 
1995) than R, the core's radius. 

Basing upon the realization that a core-crust system's temperature is generally much smaller 
than the typical chemical potentials of the constituents (u, d, s, e~), we further assume the core 
is at zero temperature. And we include first-order a c effects in our calculations. 

The thermodynamic potentials (per unit volume) as functions of the chemical potentials of 
the constituents read (Farhi & Jaffe 1984; AFO): 

r> / \ 1 r / 2 2 4\l/2/ 2 ^ 2 4\ , ^ 4 81 /Ms + (Ms 2 ~ TO^C 4 ) 1 2 

n s {fi s ) = - ^o lAMMs -m s c) 1 (p s - -m s c ) + -m s c 5 ln( ) 

2 a c ro/ / 2 2 4\l/2 2 4i /Ms + (Ms 2 — TO^C 4 ) 1 / 2 2 0/ 2 2 i\2 

[3(^ s (m s ~ rris VT /Z - m s z c^ln( )y - 2(/j,/ - m s z c^) z 

ir m s c z 

-3m s 4 c 8 ln 2 (^) + 6M— )(M s ™* 2 c 4 (Ms 2 ~ m^c 4 ) 1 / 2 
M« Ms 
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—m s c ln( ^ ))]}, (3) 

The renormalization point for the SQM, pn, which appears in Eq. @, is chosen to be 313 
MeV in this work, because of the reasons pointed out by AFO. The number densities for every 
constituent can be expressed in terms of pi (i = u,d,s,e) from (Farhi & Jaffe 1984;AFO) 

ni(Pi) = -^ 1 , (5) 



and charge neutrality requires 



^n M - -n d - -n s - n e = . (6) 



In fact, Eq. could be rendered to an equation of pi (i = u, d, s, e~) if we substitute Eqs. 
into it. Later in this section, we will demonstrate that all pi at a given energy density 
can be determined if combined with the chemical equilibrium conditions and the energy density 
equation, and we will point out that V q is equal to p e . 

Chemical equilibrium between the three quark flavors and the electrons is maintained by 
weak interactions (i.e., /3-stable SQM, Farhi & Jaffe 1984; AFO) 

d -> u + e+V e , (7) 

u + e — > d + v e , (8) 

s — > u + e + V e , (9) 

u + e — > s + v e , (10) 

and 

s + u <r+ d + u . (11) 

The loss of neutrinos by the star implies that their chemical potential is equal to zero (we 
ignore the effect brought by the finite mass of neutrinos). Hence, at equilibrium the chemical 
potentials should obey: 

H d = H s =H, (12) 

fl u + /J, e = (J, . (13) 

Combined with the condition of charge neutrality (Eq. J^J), these equations leave us with only 
one independent chemical potential, which we have denoted by [i (see Eq. H120 ). Therefore, we 
need still one equation to make the system close. 
The total energy density p is given by 

p = y"^(f2i(/Xj) + fijiiiini)) + B, i = u,d,s,e~. (14) 

Since all p,i could be expressed in terms of /i, solving Eq. (|14l) can determine p, if given the value 
of p, and so forth, pi. 

If the Thomas-Fermi model is invoked to describe the electrons associated with the quarks, 
we could easily arrive at the conclusion that V = p e , i.e., the electrostatic potential is equal to 
the chemical potential for the electrons (see, e.g., AFO; Kettner et al 1995). The derivation is 
as follows. The number density of electrons is given by the local Fermi momentum P e , 

n e = P e 3 /(3TT 2 h 3 ). (15) 
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On the other hand 



n e = -dfl e /dii e = Me 3 /(37r 2 7i 3 c 3 ). 



(16) 



From the above two equations, we get fi e — P e c. Since the electrons are confined within a 
sphere ol infinite radius, their total energy —V + Pc should obey, —V + Pc < — V(oo) — 0, i.e., 
— V + P e c = 0, and so forth, 



As fj, e at R m could be calculated following the procedure stated above, V q is readily at hand. 
3.2 Numerical Results 

Since p, e decreases with density (Kettner et al 1995), which reflects the fact that less electrons 
are needed in denser SQM, the electrostatic potential of electrons increases monotonically from 
the center toward the surface of strange stars. Because when r < R m the electrical charge 
neutrality is locally satisfied, only the electron electrostatic potential at R m , corresponding 
to V q in Eq. Q is responsible for supporting the nuclear crust. Owing to the equation of 
state P = (p — 4B)/3 (Witten 1984) 1 , the energy density at the surface of SQM cores is 
universally equal to AB, independent of the star's mass, or in other words, of the central density. 
Hence, although in our calculations we fixed the energy density at 4J3, our results will stand 
for the complete equilibrium sequences of compact SQM core-crust systems as determined by 
Glendenning, Kettner, & Weber (1995). 

There exist large uncertainties in the three QCD-related parameters: the s quark mass, m s , 
the strong interaction coupling constant, a c , and the bag constant, B. Our calculations are 
actually carried out in a space expanded by those parameters. The exact values of them are 
unknown but are probably constrained within: 50 MeV < m s < 340 MeV, < a c < 0.6, and 
135 MeV < B 1 / 4 < 165 MeV. 

The contours at fixed V q in the a c — m s plane for B 1 / 4 = 135, 145, 155, 165 MeV are 
respectively presented in Fig. l(a)-(d). The solid and dashed curves refer to V q =26 MeV (the 
electron chemical potential at which neutron drip occurs) and V q =0 MeV, respectively. The 
gray regions show where the energy per baryon of the SQM, /x n (/x n = fi u + /z<j + fi s ) , exceeds 
the lowest energy per baryon found in nuclei, which is 930 MeV for 56 Fe. The black regions are 
regions of no physical solutions. 

To support a crust at a density leading to neutron drip, V q must, at least, be larger than 
26 MeV (in order to establish an outwardly directed field, or in other words, to get meaningful 
solutions to Eq. (JTJ). So, only parameters in regions above the solid curves can make such 
crusts possible, and, of course, they should avoid the gray and black regions. Conventionally, 
the three parameters are chosen to be m s =200 MeV, a c =0.3, and B 1 ' i =145 MeV. The cross 
in Fig. 1(6) represents this set. Obviously, it lives outside the permitted region, and actually V q 
equals roughly 20 MeV. We point out emphatically that V q < 26 MeV here. 

Fig. 1 clearly shows that for QCD parameters within rather wide ranges, owing to the 
insufficient electrostatic potential of electrons, an SQM core is not capable of carrying a crust 
at pdrip suspended. 

Another factor we should take into account is the gap width, zq. Note that the lattice 
spacing in the crust is ^200 fm, in the same order of zg, while the model (i.e., Eq. (QJ) assumed 
a smooth distribution of the ionic charges in the crust. In order to make their analysis self- 
consistent, AFO assumed that the width of the gap should be larger than 200 fm. zq is entirely 
determined by V c and V q in AFO's model, we thus can also determine V q if given the values 

1 Although this EOS is derived in the limit m s — > 0, a c — » 0, for intermediate values of m s this 
equation is less than 4% different from the full expression (AFO). 



V = P e c = fx e . 



(17) 
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of zg and V c . If zq— 200 fm and V c = 26 MeV (viz. a neutron-drip crust), then V q was found 
to be ^125 MeV after solving Eq. with the associated boundary conditions (i.e., V — > 
as 2 — > — oo, V — ► 14 as z — > +oo, and V,dV/dz both continuous at z=0 and at %=zq) as 
well as the approximation dV/dz — for z > zq (viz. the gravitation forces are neglected 2 ). 
Therefore, if a gap width larger than 200 fm is really a necessary qualification for the stability 
of a core-crust system, the assumed presence of crusts at neutron drip density would have been 
completely ruled out, since no QCD parameters can make V q so high, at least within the ranges 
in which our calculations were done. 

Although our calculations were limited to SQM at zero temperature, we can state unam- 
biguously that a finite temperature can only make such dense crusts' situation even worse. 
According to Kettner et al. (1995)'s calculations (see Fig. 12 in their paper), V q decreases with 
temperature of the SQM. It thus becomes even more unlikely that a neutron-drip crust can 
exist above strange stars at higher temperatures. 

4 CONCLUSIONS 

We have studied the electrostatic potential of electrons near the surface of SQM cores in a 
parameter space expanded by m s , a c and B, in order to evaluate the feasibility of the existence 
of crusts at neutron drip density on strange stars. Our numerical results indicate: 

1. For QCD parameters within rather large ranges (Fig. 1), the electrons in an SQM core 
are incapable of establishing an outwardly directed electric field to carry a neutron-drip crust 
suspended. 

2. If it is a sound criteria that the gap width zq should be larger than 200 fm, then the 
possibility of neutron-drip crust's presence on strange stars is completed ruled out by our 
calculations. 

Therefore, our results prefer against the presence of such crusts in nature. 
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Fig. 1 Contours of fixed V q in a c - m s plane for B 1 / 4 = 135, 145, 155, 165 MeV. 
The cross near the center of (b) marks the position of conventional choice of QCD 
parameters, i.e., to s =200 MeV, a c =0.3, and B 1 ' A =1A5 MeV. The solid and dashed 
curves refer to V q =2Q MeV (the electron chemical potential at which neutron drip 
occurs) and ^=0 MeV, respectively. The gray regions show where the energy per 
baryon of the SQM exceeds the lowest energy per baryon found in nuclei, which is 930 
MeV for iron. The black regions are regions of no physical solutions. Obviously, for 
QCD parameters within rather wide ranges (including the conventional choice), SQM 
cores are incapable of supporting crusts at neutron drip density (see text for details) . 
These results are independent of stellar mass, or in other words, of central density, 
because in the MIT bag model the energy density at zero pressure (corresponding to 
the surface of SQM cores) is universally equal to AB. 



